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We study the deconfinement and chiral restoration transitions in the context of non-local 
0^ ' PNJL models, considering the impact of the presence of dynamical quarks on the scale 

, parameter appearing in the Polyakov potential. We show that the corresponding critical 

C^Hi temperatures are naturally entangled for both zero and imaginary chemical potential, in 

good agreement with lattice QCD results. We also analyze the Robcrgc Weiss transition, 

(D ■ 

' which is found to be first order at the associated endpoint. 
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I The detailed understanding of the behavior of strongly interacting matter at finite temperature 



and baryon density represents an issue of great interest in particle physics jl| . From the theoretical 
point of view, this problem can be addressed through lattice QCD calculations [^tJ], which have 
been significantly improved in the last years. However, this ab initio approach is not yet able to 
provide a full understanding of the QCD phase diagram. One well-known difficulty is given by the 
so-called sign problem, which arises when dealing with finite real chemical potentials. Thus, it is 
worth to develop alternative approaches, such as the study of effective models that show consistency 
with lattice QCD results and can be extrapolated into regions not accessible by lattice techniques. 
One of these effective theories, proposed quite recently, is the so-called Polyakov-Nambu-Jona- 
Lasinio (PNJL) model 5-11], an extension of the well-known NJL model [3] in which quarks are 



coupled to the Polyakov loop (PL), providing a common framework to study both the chiral and 



deconfinement transitions. As a further improvement over the (local) PNJL mode 



, extensions that 



13l4l5l]. The non-local 



include covariant non-local quark interactions have also been considered 
character of the interactions arises naturally in the context of several successful approaches to 
low-energy quark dynamics, and leads to a momentum dependence in the quark propagator that 
can be made consistent [3] with lattice results. It has been shown 17-20] that non-local models 
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provide a satisfactory description of hadron properties at zero temperature and density. Moreover, 



it has been found that, under certain conditions, it is possible to derive the main features of non- 



local PNJL models starting directly from QCD 
models have also been recently considered 22l |. 



Related Polyakov-Dyson-Schwinger equation 



The aim of the present work is to analyze the relation between the deconfinement and chiral 
restoration transitions at both zero and imaginary chemical potential /i in the context of non-local 
chiral quark models. One of the problems of the standard (local) PNJL model is that once the 

Q 

PL potential is adjusted to reproduce the pure gauge lattice QCD results, it is found [a] that the 
critical temperature for the chiral and deconfinement transitions at vanishing chemical potential, 
Tc ~ 220 MeV, is somewhat too high in comparison with the presently most accepted lattice 
result, namely Tc = 173(8) MeV for two light flavors A solution to this difficulty follows 
from the observation made in the context of the Polyakov-Quark- Meson model [23|], where it is 
claimed that in the presence of dynamical quarks one should decrease the parameter Tq which sets 
the scale of the PL potential. However, in contradiction to lattice results, in the PNJL model 
this sort of rescaling leads to a rather noticeable splitting between the deconfinement and chiral 
restoration t emp eratures. This splitting can be avoided only after the inclusion of extra eight-quark 
interactions 2J], or by assuming that the quark-quark coupling constant is some ad-hoc function 



of the Polyakov loop 



25l |. Here we show that in the case of the non-local SU(2) PNJL model 



the critical temperature can be made naturally compatible with lattice QCD estimates, without 
spoiling the entanglement between deconfinement and chiral restoration transition temperatures 
even for imaginary chemical potential. It should be stressed that the extension to imaginary 
chemical potential deserves significant theoretical interest, since lattice calculations j26l428l| become 



free of the sign problem and the corresponding results can be compared with effective model 
predictions. Moreover, the behavior in the region of imaginary chemical potential is expected to 
have implications on the QCD phase diagram at finite real values of u. Lattice calculations, as 



29|, suggest a close relation 



well as analyses based on the exact renormalization group equations 
between the deconfinement and chiral restoration transitions for imaginary chemical potentials. 
Thus, we extend our study of these transitions to the region of imaginary /i, where we also analyze 
the characteristics of the so-called Roberge- Weiss (RW) transition [s^, in particular, at the RW 
endpoint. 



Let us briefly describe the model under consideration, namely a non-local SU(2) chiral quark 
theory that includes couplings to a background color gauge field. The Euclidean effective action is 
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given by [31 1 
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ja{x)jaix) - jp{x)jp{x) + 



(1) 



where ^/^ is the Nf = 2 fermion doublet "0 = {u,d)'^ , rh = niq l2x2 is the current quark mass matrix 
in the isospin hmit, and = — iA^ is a covariant derivative, being color gauge fields. The 
nonlocal currents ja{x), jp{x) are given by 



ja{x) 
3p{x) 



(fz Q{z) ^[x + ^ Ta V (a; - ^ 



F{z) ^[x + - 



(2) 



where Ta = (l,i75T), and the functions J-{z) and ^(2;) are non-local form factors that characterize 
the interactions. Notice that even if we take for convenience the same coupling parameter Gs for 
both interaction terms, the relative strength between them is controlled by the mass parameter Kp. 

To proceed we perform a standard bosonization of the theory, introducing bosonic fields (Ti^2{x) 
and '/ra(x), and integrating out the quark fields. We will work within the mean field approximation 
(MFA), in which these bosonic fields are replaced by their vacuum expectation values ai^2 and 
TTfl = 0. Since we are interested in studying the characteristics of the chiral phase transition, we 
extend the bosonized effective action to finite temperature T and chemical potential ^u. This can be 
done by using the Matsubara formalism. Concerning the gauge fields A^, we assume that quarks 
move on a constant background field tj) = Ai = iAq = ig 6^q Ga^"'/2, where Ga are the SU(3) color 
gauge fields. Then the traced Polyakov loop, which is taken as order parameter of confinement, is 
given by $ = iTr exp(i(/)/r). We wih work in the so-called Polyakov gauge, in which the matrix 
(f) is given a diagonal representation (j) = (j)sXs + (/'s-^s- This leaves only two independent variables, 
(j)3 and (ps- The mean field traced Polyakov loop reads then 



* = i 

3 



exp 



2i (p8\ , „ f i 4>8\ ( 03 



(3) 



Within this framework the mean field thermodynamical potential $7^^^ at finite temperature 
and arbitrary (in general, complex) chemical potential is given by 







MFA 



c=r,g,b n=—oo 



(2^)3 

Here, M{p) and Z(p) are given by 



In 



+ 



^1 + S ^2 

2Gs 



+ ZY($,$*,T). (4) 



M{p) = Z{p) [m, + g{p)] , Z{p) = [1 - ^2 /(p)]"^ , 



(5) 
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where g{p) and f{p) are Fourier transforms of G{z) and J-{z). We have also defined 



Pn^p 



{2n + 1)-kT - ifi + • 



where the quantities 0c are given by the relation 
4'g = -03 + 4>8/V^, A = -^(ps/V^- 



diag( 



h), i-e. 



(6) 

03 + (ps/VS, 



To proceed we need to specify the explicit form of the Polyakov loop effective potential 
U{<I>,<^*,T). Following Ref. fl] we take 



,$*,r) = |- ^a(T)$$* + b{T) ln[l-6$$* + 4$3 + 4($*)3-3(«>«>*)2]| 



(7) 



(8) 



where the coefficients are parameterized as 

a{T) = ao + ai + a2 b{T) = 63 

The values of and 63 are fitted Q| to lattice QCD results, which in absence of dynamical quarks 
lead to a deconfinement temperature Tq ~ 270 MeV. However, as mentioned above, it has been 
argued [23| that in the presence of light dynamical quarks this value has to be modified accordingly, 



e.g. To ~ 208 MeV for Nj = 2 and Tq ~ 180 MeV for Nj = 3. Effects of this change in Tq will 
be discussed below. In addition, it is seen that 17^^^ turns out to be divergent, thus it has to be 
regularized. Here we use the same prescription as e.g. in Ref. 



32], namely 



(9) 



where Q.^"^^^ is obtained from Eq. by setting cri = CJ2 = 0, and ^^^^g regularized expression 

for the quark thermodynamical potential in the absence of fermion interactions: 







free 
reg 



-4T 



(2^ 



c=r,g,b s=±l 



exp 



(10) 



with Ep = yp"^ + iTT-q ■ The last term in Eq. Q is just a constant fixed by the condition that rj^g^ 
vanishes at T = = 0. 

The mean field values (T12 and 03.8 can be obtained from a set of four coupled "gap" equations 
that follow from the minimization of the regularized thermodynamical potential, 



nj^MFA 



QqMFA 



QqMFA 



. 



(11) 



dai da2 903 908 

Once the mean field values are obtained, the behavior of other relevant quantities as functions of 
the temperature and chemical potential can be determined. We concentrate in particular in the 
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chiral quark condensate (qq) = dQ,^g^/dmq, which together with the modulus of the Polyakov 
loop |<I>| will be taken as order parameters of the chiral restoration and deconfinement transitions, 
respectively. For simplicity, the associated susceptibilities will be defined as Xcond = d{qq)/dT and 
XPL = d\^\/dT. 

In order to fully specify the model under consideration we have to fix the model parameters as 
well as the form factors g{q) and f{q) that characterize the non-local interactions. Here we consider 
three parameter sets A, B and C, which have been introduced in Ref. [3]. Set A corresponds to 
the relatively simple case in which there is no wave function renormalization (WFR) of the quark 
propagator, i.e. f{p) = 0, Z{p) = 1, and g{p) has a Gaussian behavior, g{p) = exp (— p^/Aq). In 
set B we consider a more general case that includes quark WFR, taking also an exponential shape 
for the corresponding form factor, f{p) = exp (— p^/A|). Finally, in set C we take Lorentzian-like 
form factors, chosen in such a way that one can well reproduce the momentum dependence of mass 
and WFR functions obtained in lattice calculations. The parameter values for sets A, B and C, 
together with the corresponding predictions for several meson properties, can be found in Ref. [3]. 

Let us analyze the deconfinement and chiral transitions at vanishing chemical potential in the 
framework of the model presented above. Taking Tq as a parameter, and solving numerically 
Eqs. (Ilip for sets A, B and C, it is found that both the deconfinement and chiral restoration 
temperatures are coincident in a wide range of values of Tq. This is illustrated in Fig. 1, where we 
show the behavior of the relevant order parameters and the corresponding susceptibilities for the 
lattice inspired parameterization set C. We consider three characteristic values Tq = 270, 208 and 
180 MeV, corresponding to the presence of 0, 2 and 3 dynamical fermions, respectively [23(. It is 



clear that for this set both transitions are crossover-like, and they occur at basically the same critical 
temperature, as it is indicated by the peaks of the corresponding susceptibilities. One might notice 
that as long as Tq decreases the chiral susceptibility tends to become asymmetric around Tc, being 
somewhat broader on the high temperature side. Though this could be considered as an indication 
that for smaller values of Tq chiral symmetry is restored at a slightly higher temperature, even at 
Tq = 180 MeV the splitting between the main peak and what might be considered as a second 
broad peak is less than 10 MeV. In addition, it is worth to point out that although the coincidence 
of the deconfinement and chiral restoration critical temperatures holds for all three parameter sets 
A, B and C, the character of the transitions may be different from one another. This is shown in 
Fig. 2, where we plot the values of the critical temperatures as functions of Tq for sets A, B and 
C. We see that for set A, which does not include WFR, the transition becomes a first order one 
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for values of Tq below ~ 235 MeV. On the other hand, for the exponential parameterization with 
WFR, set B, this happens at a lower value Tq ~ 190 MeV. Finally, as already mentioned, for the 
lattice inspired parameterization set C the transitions are crossover-like for all considered values 
of Tq. It should be stressed that for Tq = 208 MeV (corresponding to our Nj = 2 model) the 
resulting critical temperatures are in good agreement with lattice QCD estimates. Indeed, we get 
Tc{0) = 173, 171 and 173 MeV for sets A, B and C, respectively. It is also important to remark 
that the nature of deconfinement and chiral restoration transitions for two light flavors in lattice 
QCD is still under debate. While most studies [3314371] favor a second order transition in the chiral 



39|. Given that in the context of the 



limit, there are also claims for a first order transition 
present non-local model the parameterizations that include WFR appear to be more realistic, the 
second order scenario turns out to be preferred. 

We consider now the situation at nonzero imaginary chemical potential. As it is well known, 
Roberge and Weiss found jsol that the thermodynamical potential of QCD in presence of an 
imaginary chemical potential fi = i9T is a periodic function of 6 with period 27r/3. This means 
that QCD is invariant under a combination of a Z3 transformation of the quark and gauge fields 
and a shift 9^9 + 2 kir/S, with integer k, in the chemical potential. Recently, it has been shown 
that this so-called extended Z3 transformation is also a symmetry of the local Polyakov-Nambu- 
Jona-Lasinio model [3]. Indeed, in the context of this model the thermodynamical potential is 
invariant under the transformations 



^'(6') ^ $(6')exp(-i2A:7r/3) , ^*{9) ^ ^*{9)ex.p{i2kTT/3) , 9^9 + 2k7r/3. (12) 

The RW periodicity is a remnant of the Z3 symmetry in the pure gauge limit. In QCD with 
dynamical quarks, if the temperature is larger than a certain value Tfjiy it can be seen that three 
vacua appear. These vacua can be classified through their Polyakov loop phases, given by ip, 
if + 2^/3 and tp + 47r/3. Roberge and Weiss showed that for T > Trvk there is a first order phase 
transition at 9 = ir/S mod 27r/3, in which the vacuum jumps to one of its Z^ images. This is known 
as the "Roberge- Weiss transition" , and the point at the end of the RW transition line in the (T, 9) 
plane, i.e. {T,9) = (Trvi/j 2vr/3), is known as the "RW endpoint". The order of the RW transition 
at the RW endpoint has been subject of considerable interest recently in the framework of lattice 
QCD 
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- 43( 1 due to the implications it might have on the QCD phase diagram a finite real fx. 
According to lattice calculations, it appears that for two light flavors the RW endpoint is first 
order for realistically small values of the current quark mass. Following these considerations it is 
important to check whether the thermodynamical potential of the non-local PNJL models studied 
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in this work does respect the extended Z3 symmetry. In fact, it is easy to show that this is the 
case. The last two terms in Eq. ^ are obviously invariant under the transformations in Eq. (I12p . 
whereas to check the invariance of the first term it is convenient to write these transformations in 
the equivalent way 

MO)^MO), M0)^M0)-'^k7rT/V3 , e^e + 2k7r/3. (13) 

Thus it can be easily proven that any sum of the form X]c=r g 6 -^[(Pn p)^]' where F is an arbitrary 
function, turns out to be invariant under the extended Z3 transformations. The invariance of the 
terms introduced in the regularization procedure [c.f. Eq.([9])] can be shown in the same way. 

Having checked that our non-local PNJL models possess the extended Z3 invariance we turn 
now to the results of the numerical analysis of the behavior of the different order parameters as 
functions of T and 9 using the value Tq = 208 MeV corresponding to two light flavors. We first 
keep T fixed, verifying that the expected periodicity of the different thermodynamical quantities 
as functions of 9 is indeed satisfied. Moreover, for T > we find the mentioned RW first order 
phase transition at 9 = tt/S, which is signalled by a discontinuity in the phase of the Polyakov 
loop field. The values obtained for Trw are 191 MeV, 188 MeV and 191 MeV for sets A, B and C, 



respectively, in good agreement with the lattice QCD estimate Tr\y = 185(9) [28|]. Concentrating 
on the sector < ^ < vr/S we observe that for values of the temperature Tc(0 = 0) <T< Tuy/ the 
order parameters for both deconfinement and chiral symmetry show signals of a phase transition at 
a given value of 9. This is clearly seen in Fig. [3] where we plot the behavior of order parameters and 
susceptibilities as functions of T taking now 9 fixed at two representative values = vr/G and 7r/3. 
The plots correspond to parameter set C. We note that while for = vr/G both deconfinement and 
chiral restoration are crossover-like, they are first order for 9 = 7r/3. As in the case of = (see the 
curves corresponding to Tq = 208 MeV in Fig. [T]), for both values of 9 the deconfinement and chiral 
restoration transitions occur at the same temperature, given by the peaks of the susceptibilities or 
the positions of the discontinuities. Although one might argue that there is a certain tendency of 
the chiral susceptibility to decay more slowly or even, in the case oi9 = vr/S, to display a very broad 
peak on the high temperature side, one can hardly conclude that both transitions apart from each 
other more than about 20 MeV, even for 9 = 7r/3. The situation is quite similar for parameter set 
B (which also includes WFR), whereas for set A one gets just first order transitions for all values 
of 9 in the range of temperatures considered. The dependence of the critical temperature Tc as a 
function of 9 is shown in Fig. U] for our three parameter sets. For comparison, we also show the 
corresponding lattice results given in Ref. [28|, which include an error of about 10% due to the 
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uncertainty in the lattice determination of Tc{0 = 0). As already mentioned, while for set A both 
the deconfinement and chiral restoration transitions are always first order, in the case of sets B 
and C there is a critical value 9 ~ 0.7 x vr/S below which the transitions become crossover- like. 
Thus, we find that for all three parameterizations the corresponding transition lines are first order 
when they meet the RW endpoint. This implies that the RW endpoint is a triple point, the RW 
transition being also first order there. The character of the RW transition at the RW endpoint is 
clearly seen in Fig. [5l where we plot the behavior of the phase of the extended Polyakov loop ^ at 
9 = 7r/3 as a function of T (the figure corresponds to set C). The extended Polyakov loop, defined 
by ^' = exp(z0)$, is by construction invariant under the transformations in Eq. (I12p . and its phase 
%[) can be taken as order parameter of the RW transition [4^. 

In summary, we have considered the impact of the feedback of the dynamical quarks on the 
parameter Tq appearing in the Polyakov potential, as proposed in Refs. [231 ]. This has been done 
here in the context of non-local PNJL models, considering three different types of non-local form 
factors. We have studied the deconfinement and chiral restoration transitions, determining the 
corresponding critical temperatures and the character of the transitions. The results are found to 
be in agreement with those obtained in lattice QCD, showing a natural entanglement between both 
critical temperatures for zero and imaginary chemical potential. We have also analyzed the Roberge 
Weiss transition, which is found to be first order at the RW endpoint for all three parameter sets. 

This work has been partially supported by CONICET (Argentina) under grants jj^ PIP 00682 
and PIP 02495, and by ANPCyT (Argentina) under grant # PICT07 03-00818. 
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FIG. 1: Order parameters for the deconfinement and chiral restoration transitions (upper panel) and cor- 
responding susceptibilities (lower panel) as functions of T for three characteristic values of To, namely 
To = 270, 208 and 180 MeV. Curves correspond to Set C. 
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FIG. 2: Critical temperature as a function of Tq for our parameter sets A, B and C. Solid and dashed lines 
stand for first order and crossover-like transitions, respectively. Note that below To ~ 190 MeV the curves 
for sets B and C almost overlap. However, while for set B the transition in that region is first order (solid 
line), for set C it is crossover-like (dashed line). 
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FIG. 3: Order parameters for the deconfinement and chiral restoration transitions (upper panel) and corre- 
sponding susceptibilities (lower panel) as functions of T ioi = tt/Q and 9 ~ 7r/3. Curves correspond to Set 
C. 
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FIG. 4: Critical temperature as a function of 9 for parameter sets A (left), B (center) and C (right). Solid 
and dashed lines stand for first order and crossover-like transitions, respectively. Dots correspond to lattice 
QCD results. Vertical solid lines correspond to the first order RW transition. 
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FIG. 5: Phase of the extended Polyakov loop ^' = exp(i0)$ as a function of T for = 7r/3. Curves 
correspond to set C. 



